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Abstract. In a fractional Cauchy problem, the usual first order time derivative is 
replaced by a fractional derivative. The fractional derivative models time delays in 
a diffusion process. The order of the fractional derivative can be distributed over 
the unit interval, to model a mixture of delay sources. In this paper, we provide 
explicit strong solutions and stochastic analogues for distributed-order fractional 
Cauchy problems on bounded domains with Dirichlct boundary conditions. Sto- 
chastic solutions are constructed using a non-Markovian time change of a killed 
Markov process generated by a uniformly elliptic second order space derivative op- 
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1. Introduction 



This paper develops explicit strong solutions for distributed-order fractional Cauchy 
problems on bounded domains D C M. d with Dirichlet boundary conditions. Cauchy 
problems du/dt = Lu model diffusion processes. The simplest case L = A = 
">") • d 2 u/dx 2 governs a Brownian motion B(t) with density u(t, x), for which the square 

root scaling u(t,x) = t~ l l 2 u{l,t~ l l 2 x) pertains [19] . The fractional Cauchy problem 
d^u/df 13 = Lu with < (3 < 1 models anomalous sub-diffusion, in which a cloud of 
particles spreads slower than the square root of time [23 ESI E3 SSI HZ] • When L = A, 
the solution u(t,x) is the density of a time-changed Brownian motion B(E t ), where 
the non-Markovian time change E t = inf{r > : W T > t} is the inverse, or first pas- 
sage time, of a stable subordinator W t with index (5. The scaling W ct = c l ^W t in law 
implies E ct = c@E t in law for the inverse process, so that u(t,x) = t"^l 2 u{\,t~^l 2 x). 
The process B(E t ) is the long-time scaling limit of a random walk "5U1 EH] when the 
random waiting times between jumps belong to the /3-stable domain of attraction. 
Roughly speaking, a power-law distribution of waiting times leads to a fractional time 
derivative in the governing equation, and its power law index equals the order of the 
fractional time derivative. For a uniformly elliptic operator L on a bounded domain 
D C M. d , under suitable technical conditions and assuming Dirichlet boundary con- 
ditions, the Cauchy problem governs a Markov process X(t) killed at the boundary. 
Then the fractional Cauchy problem governs the time-changed process X(E t ) [33] . 
Recently, Barlow and Cerny [6] obtained B(E t ) as the scaling limit of a random walk 
in a random environment, when the transition rates are drawn from a power-law 
distribution with index (3. 

Fractional derivatives are almost as old as their integer-order cousins "5(3 H2]. Frac- 
tional diffusion equations are important in physics, finance, hydrology, and many 
other areas [2TI l"2"6"[ 155*1 135] . Nigmatullin [5j5] gave a physical derivation of the fractional 
Cauchy problem. The mathematical study of fractional Cauchy problems was initi- 
ated by Schneider and Wyss [35] and Kochubei "25J [26] . Fractional Cauchy problems 
were used by Zaslavsky [J7] as a model for Hamiltonian chaos. Stochastic solutions 
of fractional Cauchy problems are the basis for particle tracking schemes [TTT [1"5| [39] . 

In some applications, the waiting times between particle jumps evolve according to 
a more complicated process that cannot be adequately described by a single power law. 
Then, a waiting time model that is conditionally power law leads to a distributed- 
order fractional derivative in time, defined by integrating the fractional derivative 
of order (3 against the probability distribution of the power- law index [52]. The 
resulting distributed-order fractional Cauchy problem provides a more flexible model 
for anomalous sub-diffusion. The Levy measure of a stable subordinator with index 
f3 is integrated against the power law index distribution to define a subordinator 
W t . Its inverse E t produces a stochastic solution X(E t ) of the distributed-order 
fractional Cauchy problem on Mr, when X(t) solves the original Cauchy problem 
du/dt = Lu. Mild solutions and stochastic solutions for distributed-order fractional 
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Cauchy problems on M. d were developed by Meerschaert and Schemer [32] , Kovacs and 
Meerschaert [2B], and Hahn, Kobayashi and Umarov [22]. Kochubei [2Z] obtained 
strong solutions for the case L = A, and proved the uniqueness for more general 
uniformly elliptic operators L. 

In this paper, we extend the basic approach of [301 EI] to solve distributed-order 
fractional Cauchy problems with Dirichlet boundary conditions on bounded domains 
D C M. d . Our recent paper [31] treated fractional Cauchy problems on bounded 
domains using eigenfunction expansions and killed Markov processes. Here, we extend 
that theory to the distributed-order fractional Cauchy problems, using the inverse 
subordinators from [2BJ E2], along with some deep technical results from Kochubei 
[21]. We construct unique classical solutions, identify the stochastic process governed 
by the distributed order fractional Cauchy problem, and prove existence. At the 
end of this paper, we discuss also some open problems in the literature, including 
extensions to jump processes. 

2. Heat kernels on bounded domains 

Let D be a bounded domain in R d . We denote by C k (D), C k ' a (D), C k (D) the space 
of k-times differentiable functions in D, the space of fc-times differential functions with 
k-th derivative is Holder continuous of index a, and the space of functions that have 
all the derivatives up to order k extendable continuously up to the boundary dD of 
D, respectively. We refer to [M] for a detailed discussion of these spaces and concepts 
in this section. 

A uniformly elliptic operator in divergence form is defined on C 2 functions by 

= V d ( a ij( X )( du / dx i)) 

with ciij(x) = ciji(x) and, for some A > 0, 



(2.2) \j2y* < E "vWwj ^ A_1 E y y G Rd - 

i=l *>i=l i=l 

The operator L acts on the Hilbert space L 2 (D). We define the initial domain 
C£°(D) of the operator as follows. We say that / is in C£°(D) if / e C°°(D) and 
f(x) = for all x G dD. This condition incorporates the notion of Dirichlet boundary 
conditions. 

If X t is a solution to 

dX t = a(X t )dB t + b(X t )dt, X = x , 

where a is a d x d matrix and B t is a Brownian motion, then X t is associated with 
the operator L with a = aa T (see Chapters 1 and 5 of Bass [7]). Define the first exit 
time as td(X) = inf{i > : X t D}. Then the semigroup defined by 

(2.3) T D (t)f(x) = E x [f(X t )I(r D (X)) > t)} 
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has generator L, which follows by an application of the ltd formula. 

Let D be a bounded domain in ~R d and L be a uniformly elliptic operator of diver- 
gence form with Dirichlet boundary conditions on D. Then there exists a constant A 
such that for all x G D, 

d 

(2.4) Z)M*)I<A, 

see, for example, Davies [T7J Chapter 6] 

Let Tz>(t) be the corresponding semigroup. Then To(i) is an ultracontractive semi- 
group (even an intrinsically ultracontractive semigroup), see Corollary 3.2.8, Theorem 
2.1.4, Theorem 4.2.4, and Note 4.6.10 in [TB]. Every ultracontractive semigroup has 
a kernel for the killed semigroup on a bounded domain which can be represented 
as a series expansion of the eigenvalues and the eigenfunctions of L (cf. [HJ The- 
orems 2.1.4 and 2.3.6] and [20| Theorems 8.37 and 8.38] ). There exist eigenvalues 
< i^i < < 1^3 ■ ■ ■ , such that \i n — > oo, as n — > oo, with the corresponding 
complete orthonormal set (in L 2 ) of eigenfunctions ifj n of the operator L satisfying 



(2.5) Ltf) n (x) = -fx n tp n (x), x E D : ^ n \dD = 0. 

In this case, 

oo 

p D (t,x,y) =^e^Vn(s)^n(2/) 

n=l 

is the heat kernel of the killed semigroup Tp. The series converges absolutely and 
uniformly on [to? oo) x D x D for all t > 0. 

Denote the Laplace transform t —> s oi u(t,x) by 

POO 

u(s,x) — / e~ st u(t, x)dt. 



Since we are working on a bounded domain, the Fourier transform methods in [30] are 
not useful. Instead, we will employ Hilbert space methods. Hence, given a complete 
orthonormal basis {vjj n (x)} on L 2 (D), we will call 



u(t,n) = / ip n (x)u(t,x)dx, 
Jd 

and 

u(s,n) = / ip n (x) / e~ st u(t, x)dtdx 



jd Jo 

(2.6) = / ijj n (x)u(s, x)dx 

Jd 

"OO 

e~ st u(t,x)dt (when Fubini's condition holds) 

o 



respectively the ip n and the ^-Laplace transforms. Since {ip n } is a complete or- 
thonormal basis for L 2 (D), we can invert the ^-transform to obtain 

u(t,x) = ^2u(t,n)ip n (x) 

n 

for any t > 0, where the above series converges in the L 2 sense (e.g., see (JU Propo- 
sition 10.8.27]). 

Suppose D satisfies the uniform exterior cone condition. Let {X t } be a Markov 
process in ~R d with generator L, and / be continuous on D. Then the semigroup 

(2.7) 

T D (t)f(x) = E x [f(X t )I( TD (X) >t)}= / p D (t,x,y)f(y)dy = Ve^V^/H 

jD n=l 



solves the Dirichlet initial-boundary value problem in D: 

= L D u(t,x), xeD,t>0, 

u(t, x) = 0, x G (9-D, 
u(0, x) = f(x), x G D. 
See [7J Theorem 6.3, page 177] or [TBi Theorem 2.1.4]. 

Remark 2.1. Let L°°(D) = {/ : WfW^ < oo}, where WfW^ = ess sup \f\. The 
eigenfunctions belong to L°°(D) fl C a (D) for some a > 0, by [2"0"| Theorems 8.15 and 
8.24]. If D satisfies the uniform exterior cone condition, then all the eigenfunctions 
belong to C a (D) by [HI Theorem 8.29]. If ay G C a (L>) and <9£> G C 1 '", then the 
eigenfunctions belong to C 1)0l (D) by [201 Corollary 8.36]. If G C°°(D), then each 
eigenfunction of L is in C°°(D) by [201 Corollary 8.11]. If ay G C°°{D) and &D G 
then each eigenfunction of L is in C°°(D) by [201 Theorem 8.13]. 

In the case L = A, the Laplacian, the corresponding Markov process is a Brownian 
motion. We denote the eigenvalues and the eigenfunctions of A on D, with Dirichlet 
boundary conditions, by {A n , n }^=i, where (p n G C°°(D). 

3. Distributed order fractional derivatives 

Fractional derivatives in time are useful for physical models that involve sticking or 
trapping [30]. They are closely connected to random walk models with long waiting 
times between particle jumps [31]. The fractional derivatives are essentially convo- 
lutions with a power law. Various forms of the fractional derivative can be defined, 
depending on the domain of the power law kernel, and the way boundary points are 
handled [3H1I12]- The Caputo fractional derivative [TU] is defined for < (3 < 1 as 

d^u{t,x) 1 f* du(r,x) dr 

(3.1) 



dtP T(1-P)J dr (t-r)l 3 ' 



Its Laplace transform 

(3.2) e~ st 9Pu ^ X) ds = S Pu(s, x) - s^uiO, x) 

incorporates the initial value in the same way as the first derivative. The Caputo 
derivative is useful for solving differential equations that involve a fractional time 
derivative [211 SO] , because it naturally incorporates initial values. 

Fractional time derivatives emerge in anomalous diffusion models, when particles 
wait a long time between jumps. In the standard model, called a continuous time 
random walk (CTRW), a particle waits a random time J n > and then takes a step 
of random size Y n . For the purposes of this paper, we may assume that the two 
sequences of i.i.d. random variables (J n ) and (Y n ) are independent. This is called 
an uncoupled CTRW. The particle arrives at location X(n) = Yi + • • • + Y n at time 
T(n) = J\ + • • • + J n . Since N t = max{n > : Tin) < t] is the number of jumps 
by time t > 0, the particle location at time t is X(N t ). If EY n = and < oo 

then, as the time scale c — > oo, the random walk of particle jumps has a scaling limit 
c~ l / 2 X(\ct\) B(t), a standard Brownian motion. If P(J n > t) ~ ct~^ for some 
< (5 < 1 and c > 0, then the scaling limit c~ l ^T{\ct\) ^ W t is a strictly increasing 
stable Levy process with index (3, sometimes called a stable subordinator. The jump 
times T{n) and the number of jumps N t are inverses {N t > n} = {T(n) < t}, and 
it follows that the scaling limits are also inverses [3H Theorem 3.2]: c~^N ct E h 
where 

(3.3) Et = inf{r : W T > t}, 

so that {E t < t} = {W T > t}. A continuous mapping argument [31], Theo- 
rem 4.2] yields the CTRW scaling limit: Heuristically, since N ct ~ c^Et, we have 
c~ p/2 X(N[ ct }) w (c^y 1/2 X(c^E t ) w B(E t ), a time-changed Brownian motion. The 
density u(t,x) of the process B(E t ) solves a fractional Cauchy problem 

d u(t,x) d 2 u(t,x) 
dt 13 dx 2 

for some a > 0, where the order of the fractional derivative equals the index of 
the stable subordinator. Roughly speaking, if the probability of waiting longer than 
time t > between jumps falls off like t _/3 , then the limiting particle density solves 
a diffusion equation that involves a fractional time derivative of the same order f3. 
Hence, the fractional derivatives in time model sticking or trapping of particles for 
long (power-law distributed) periods of time. 

A more flexible model for diffusion processes can be obtained by considering a 
sequence of CTRW. At each scale c > 0, we are given i.i.d. waiting times ( J^) and i.i.d. 
jumps (Yn)- Assume the waiting times and jumps form triangular arrays whose row 
sums converge in distribution. Letting X c (n) = Yf+- ■ -+Y° and T c (n) = Jf+- • -+</^, 
we require that X c (cu) =^ A(t) and T c {cu) =^ W t as c — > oo, where the limits A(t) 
and Wt are independent Levy processes. Letting N£ = maxjn > : T c {n) < £}, the 
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CTRW scaling limit X c (Nf) A(E t ) [331 Theorem 2.1]. A power-law mixture model 
for waiting times was proposed in [32]: Take an i.i.d. sequence of mixing variables 
(Bi) with < Bi < 1 and assume P{Jf > u\Bi = (3} = c~ x u~^ for u > c~ 1//3 , 
so that the waiting times are power laws conditional on the mixing variables. The 
waiting time process T c (cu) =>- Wt a nondecreasing Levy process, or subordinator, 
with E[e~ sM/ '] = e~ t ^ w ( a ' and Laplace exponent 



(3.4) ^ w (s) = / (e~ sx - l)4> w {dx). 

Jo 

The Levy measure 

(3.5) (f>w(t, oo) = [ r p fi(dP), 







where /i is the distribution of the mixing variable [32j Theorem 3.4 and Remark 5.1]. 
A computation [321 Eq. (3.18)] using / °°(1 - e'^fit-^dt = T(l - (3)s p shows that 

(3.6) M*)= / s^T(l-l3)fi(dl3). 

Jo 

Then c _1 iV f c => E t , the inverse subordinator [321 Theorem 3.10]. The general infinitely 
divisible Levy process limit A(t) forms a strongly continuous convolution semigroup 
with generator L (e.g., see [3]) and the corresponding CTRW scaling limit A(E t ) is 
the stochastic solution to the distributed order-fractional Cauchy problem [321 Eq. 
(5.12)]defined by 

(3.7) B {u) u(t,x) = Lu(t,x), 
where the distributed order fractional derivative 

(3.8) B^u(t, x) := C 8Pu !f' X K {df3), u{d(3) = T(l - p)n(dp). 

Jo OtP 

To ensure that is well-defined, we impose the condition 

f 1 1 

(3.9) / -^{d(3)<oo 



1-/3 

as in [3"2"l Eq. (3.3)]. Since T(x) ~ 1/x, as x — > 0+, this ensures that u(df3) is a finite 
measure on (0, 1). 

Using triangular array limits for CTRW allows a more flexible limit model. For 
example, suppose Y{ are i.i.d. with mean \i and variance a 2 , and let Y? = c~ l [i + 
c -1 / 2 (Tj — jx) so that X c (cu) A{t) a Brownian motion with drift. Then the density 
u(t,x) of the CTRW scaling limit A(E t ) solves (J33D with 

d d 2 
L = —v— — h a 



dx dx 2 
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for some a > 0. A triangular array of jumps with two spatial scales, one for the mean 
jump and another for the deviation from the mean, is necessary to get a drift term 
in the limit. 

Since 0vf(O, oo) = oo in (13. 4p . Theorem 3.1 in [33] implies that the inverse subor- 
dinator E t has density 



This same condition ensures also that E t is almost surely continuous, since Wt jumps 
in every interval, and hence is strictly increasing. Further, it follows from the defini- 
tion (13. 3p that E t is monotone nondecreasing. 

The following lemma shows that h(t, A) = E[e _A£;t ] is an eigenfunction of the 
distributed-order fractional derivative D^. 

Lemma 3.1. For any X > 0, h(t, A) = / °° e~ Xx g(t,x) dx = K[e~ XEt ] satisfies 



Proof. First note that h(0, A) = E(l) = 1. Using (13.21) and (13.81) . compute the Laplace 
transform of B^h(t,X) as 




(3.10) 



OM/i(t, A) = -Xh(t, A); h(0, A) = 1. 



(3.11) 




by applying a Fubini argument which holds because ipw(s) < °°- 
The Laplace transform of g(t,x) is given by [331 Eq. (3.13)]: 



(3.12) 
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Then the double Laplace transform 

roo roc / roc 

h(s,\):= / e~ st h(t,X)dt = / e- xt ( / e~ Xx g(t, x)dx ) dt 



oo 

e ~ XX (j e- st g(t,x)dt] 



(3.13) 
(3.14) 

That is, h(s, A) satisfies 



s JO 



s(X + ipw(s)) 



(3.15) Xh(s, A) = (- - fc(s, A))^ H/ (s). 

Since has continuous paths, the dominated convergence theorem implies that t — * 
E[ e -A£(t)] = fr(t, A) is a continuous function. Then fl3~T0|) follows from (133T1) . fl3~T5l) 
and the uniqueness of the Laplace transform. □ 

Lemma 3.2. Suppose that n(d[3) = p{(3)d(5, the function [3 \— > T(l — (3)p(f3) is in 
C^O, 1], supp(fi) = [Po,Pi] C (0, 1) and fi(Pi) ^ 0. Suppose also that 

rPi 

(3.16) C((3 0l (3 u p) = / sin(/37r)r(l-/3)p(/3)d/3>0. 
Then \d t h(t,X)\ < Xk(i), where 

(3.17) k(t) = [C , (/5 ,A,p)7r]- 1 [r(l-/3 1 )t /3l - 1 + r(l-/5o)t /3 °- 1 ]. 
In this case, h(t,X) is a classical solution to (13.101) . 

Proof. Using (2.19) in Kochubei [27], which follows from inverting the Laplace trans- 
form in (13. 13[) of h(t, A), we have 

—A r°° 

(3.18) h(t,X) = — / r _1 e-* r $(r, l)dr 

n Jo 

where 

/ 1 r^sin(/3vr)r(l-/3M/5)^ 



$ r, 1 



[Jo 1 r/ 3 cos(/57r)r(l - p)p(f3)d{3 + A] 2 + r/ 3 sin(/?7r)r(l - P)p(f3)d(3} 2 ' 



First we show that \d t h(t, A)| < oo. Note that 



(3.19) 

\d t h(X,t)\ 
X rco 



7T 



r^[d t e- tr ]§{r, l)dr 



k JO 

X rco 



e $(r, l)dr 



TV 



e~ tr Jo r 13 sin(/37r)r(l - f3)p(P)d/3 



o [/„ rP cos(/57r)r(l - p)p(p)dp + A] 2 + r/ 3 sin(/ftr)r(l - f3)p{(3)d{3\ 
e~ tr dr 



-dr 



o / Vsin(/?7r)r(l-/?)p(/?)d/3 
XI (t), (say), 



where Z(t) is a function of £ only. In the case of a simple fractional derivative this l(t) 
is given by Ct 13 ^ 1 . 
Now, 



-1 rPi 

(3.20) / sm((37r)T(l-/3)p(P)d(3> sin(/?7r)r(l - (3)p((3)d(3 = (?(#,, ft, p) > 0, 

i/3o 



by assumption (13.161) . 

For r > 1, and /3q < j3 < fli < 1, we have < r' 3 < r' 31 and so 



(3.21) 



/9i 



srin(ftr)r(l - (3)p{(3)d(3 > / sin(ftr)r(l - (3)p{(5)d(5 



r^C(A>,ft,p). 



For < r < 1, and /3q < < ft < 1, we have r"° >r@> r^ 1 and 



so 



sin(ftr)r(l - f3)p((3)d(3 > / r^ 1 sin(ftr)r(l - (3)p{(3)d(5 

(3.22) J/30 J/30 



r*C(A>,ft,p). 
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Using the above facts, we obtain 

w 



IT 



7T 



<7T- 1 



e tr dr 


/ Vsm(/37r)r(l - 




e~ tr dr 




/ Vsin(/37r)r(l - 




e~ tr dr 





e dr 



o sm{f3it)T{l - f3)p{(3)df3 Ji J£ r/»fdn(/?7r)r(l - P)p(PW 



e tr dr 



<[C{^P u p)ix}- 1 



r- pi e~ tr dr+ / r^e^dr 
Ji 



< [C(/3 ,/3i,p)7r]- 1 [r(l-/3 1 )#- 1 + r(l-/3oK°- 1 ] =k(t) 
and so \d t h(t,X)\ < Xk(t). Hence, it follows from (13. 8p that 



|D^/i(t,A)| < 



< 



< 



j tLh(t,\)r(i-p) P (p)dp 

1 dh{s,X) 



o r(l-/3)J 
i 1 



<9s 



-A- 



r(l-/3)p(/3)d/3 



ds 



T{l-(3)p{(3)dp 



/ r(i-/3)7 -^(t- a )/» 

< oo, 

using (I3.9I) and the beta density formula. Thus, the distributed-order derivative 
D^/i(t, A) is defined in the classical sense. □ 



4. Distributed order fractional Cauchy problems 

Fractional Cauchy problems replace the usual first-order time derivative with its 
fractional analogue. In this section, we prove classical (strong) solutions to distributed- 



order fractional Cauchy problems 



u 



Lu on bounded domains Del We give 



also an explicit solution formula, based on the solution of the corresponding Cauchy 
problem. Our methods are inspired by the approach in [30], where Laplace transforms 
are used to handle the fractional time derivative, and spatial derivative operators (or 
more generally, pseudo-differential operators) are treated with Fourier transforms. In 
the present paper, we use an eigenfunction expansion in place of Fourier transforms, 
since we are operating on a bounded domain. Our first result, Theorem I4.1[ is fo- 
cused on a distributed-order fractional diffusion with L = A, and we lay out all the 
details of the argument in the most familiar setting. Then, in Theorem I4.6I we use 
separation of variables to extend this approach to uniformly elliptic generators L. In 
the process, we explicate the stochastic solutions in terms of killed Markov processes. 

n 



Let -Dqo = (0, oo) x D and define 

H A poo) = R : A« e C^), 

IM*.«)I < g£L°°(D), t>o}, 

where k(t) is denned by (I3.17p . 

We will write u G C k (D) to mean that for each fixed t > 0, u(t, •) G C k (D), and 
u G C^(D oo) to mean that it G C k (D 00 ) and is bounded. Let Td(X) = inf{t > : 
X(t) (fi D} denote the first exit time of the stochastic process X = {X(t)}. 

Theorem 4.1. Let D be a bounded domain with dD G C 1,a for some < a < 1, 
and To{t) be the killed semigroup of Brownian motion {X(t)} on D. Let E t be the 
inverse (13.31) of the subordinatorW t , independent of {X (t)} , with Levy measure (13.51) . 
Suppose that fi(dj3) = p{[3)d(3, as in Lemma Iff.H and W^> is the distributed- order 
fractional derivative defined by (13.81) . Then, for any f G D(Ap) C\C 1 (D) C\C 2 (D) for 
which the eigenfunction expansion ( of Af ) with respect to the complete orthonormal 
basis {(j) n : n G N} converges uniformly and absolutely, the unique (classical) solution 
of the distributed order fractional Cauchy problem 

(4.1) B {u) u{t,x) = Au(t,x); x G D, t > 

u(t,x) = 0, x G dD, t>0, 
u(0,x) = f(x), x G D, 

for u G 7Ya(Ax>) n C 6 (Ax>) n C l (D), is given by 

u(t,x) = E x [f(X(E t ))I(r D (X)>E t )} 

POO 

(4.2) = / T D (l)f(x)g(t,l)dl 

Jo 

oo 

(4.3) = ^/(n)0 n (x)/i(t,A n ), 

n=l 

where h(t,\) = K(e~ XEt ) = / °° e~ Xx g(t, x) dx is the Laplace transform of E t . 
Proof. Assume that u(t,x) solves (14. ip . Using Green's second identity, we obtain 

d(t>r, , du 



uA(p n — 4> n Au\dx = 

D J 3D 



li- 



ds = 0, 



36 ^' l d9 

since u\qd = = (f) n \dD an d u, <p n G C 1 (/)). Hence, the n -transform of Am is 

4> n {x)Au{t, x)dx = I u{t, x)A(j) n (x)dx = — X n / u(t,x)(p n (x)dx = —X n u(t,n), 

D J D J D 

as cj) n is the eigenfunction of the Laplacian corresponding to eigenvalue X n . 
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Next we need to show that the (j) n transform commutes with DM. We need to show 
that we can interchange derivatives and integrals as follows. Observe that 

<f) n {x)B {u) u{t,x)dx 



i) 



0nO) / Q^u{t,x)r(l- P)p(P)d/3dx 



4>n(x) 



D 



du(s,x) ds 



(4.4) 



o T(1-(3)J ds (t-sy 
'* du(s, x) ds 



r(l - p)p(p)dpdx 



D 




J 



1 ft 




ds (t - s)P 







p((3)d(3dx 



'0 Jo \jd 

»i ft 



<Pn{x)—u{s, x)dx 



ds 




JO 

1 



d_ 

ds 
1 



4> n (x)u(s, x)dx 



D 



(t - s y 

ds 



(t - 8 y 



p(/3)dj3 (by Fubini, see below) 
p(P)dP 



T(1-P)J ds 
( u 'u(s, n). 



1 d ds 

u(s, n) . * r(l - f3)p{P)d/3 



(t - s y 



The Fubini- Tonelli argument for the interchange of order of integration in (14.41) can 
be justified as follows: 



D 



(j) n (x)I} iu) u(t J x)dx 



'* du(s, x) ds 



D 




JO 

1 ft 



ds (t - s y 



< / \<f>n{x)\ 



I) 




'0 Jo 

< I \4> n (x)\\g(x)\dx 



du(s, x) 


ds 


ds 


(t - 8 y 



1 ft 



D 




p{f3)d(3dx 
p((3)dpdx 
p(j3)d/3 



< V\D\m\L H D)\\g\h 



ff 

'0 Jo 

x 



[C(p , A, P )7r]- 1 [r(i - fa)^- 1 + r(i - a,)**- 1 ] 

ds 



(t - s y 

using (13.171) . Further, using the property of beta density, for < 7,77 < 1, 



(t-sp 



^ds = t^ / (l-u)^- 1 u"- 1 dM = S(l-7,77)^- 7 , 
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where B(a, b) denotes the usual beta function. Thus, 



<p n (x)B { " ) u(t,x)dx 



' D 



< V / ^D|||0n||L2(D)||fi'||L°°[C , (^o, Pl,P)ir] 1 



1 rt dt 



JO 



(t - s y 



p{(5)d(5 



+ 



1 r t 



r(i-A))a 



/5o-l 



ds 



-M(3)d(3 



{t-sy 

= \/^D|||0n.||L2(D)||fi'||L-[C(/3o,/5l,P)vr] _1 

+ r(l-/3 ) / t*>-f>B(l-l3, Po)p(J3f)dl3 
Jo 

< oo, 



which justifies the use of Fubini-Tonelli theorem in (I4.4p . 
Thus, applying the n -transforms to (14. ip . we get 



r(l-ft) / &-f>B(l-W x )p{0)dp 



(4.5) 



DM 



u(t,n) = —\ n u(t,n). 



Since w is uniformly continuous on C([0, e]) x D), it is uniformly bounded on [0, e] x D. 
Thus, by the dominated convergence theorem, we have lim^ f D u(t,x)<ft n (x)dx = 
f(n). Hence, u(0,n) = f(n). A similar argument shows that t i— > u(t,n) is a continu- 
ous function of i G [0, oo) for every n. Then, taking Laplace transforms on both sides 
of fT4~5D . we get 



(4.6) 
which leads to 
(4.7) 



[s p u(s, n) - s^ 1 u(0, n))r(l - /3)p(j3)d/3 = -X n u{s, n) 



u(s, n) 



f(n)J ( ;s^ 1 T(l-P)p(P)dP 



Recall Jq 1 s^T(l - (3)p((3)d(3 = ifay(s). Then, from 

f(n)ifj w (s) 



u[s,n} 



S(i>w(s) + A n ) 

-f(n)4>w(s) / e -^ w{s)+Xn)l dl 



Jo 

OO -| 

e" A?l 7(n)-?Ms)e-^ (s) d/, 



o 



using the property of the exponential density. 
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The (^-transform of the killed semigroup Tjj(l)f(x) = Ylm=i e ~ Xml( t ) m{x)f{m) from 
(12.71) is found as follows. Since {<f> n , n e N} is a complete orthonormal basis of L 2 (D), 
we get 



[Tn(l)f}( 



n 



D 



D 



D 



(4.9) 



<p n (x)T D (T)f(x)dx 
(j>n{x) [ p D (l,x,y)f(y)dydx 

oo 

e~ Xml <j) m (x)<j) m (y)f(y)dydx 

(y)f(y)dydx 



x 



D 



m=l 



D 



e Am ^ m (x) 



D 



m=l 

» OO 

/ <j> n (x) E e " Am Vm(x)/(m)dx 

m=l 

OO 



(j) n (x)4>m(x)dx 



m=l 
-IX 



= e- lx -f{n). 

Since To(t) is a contraction semigroup on L 2 (D), T D (t)f e L 2 (D) and hence Fubini- 
Tonelli applies. 
By (3.20) in [32], we have 



-^(s)e"^ (s)/ = / e- st g(t,l)dt 



(4.10) 



where g(t, /) is the smooth density of E t . 
Using the results (TCHD, fl4TTU|) and (USD, we get 



e~ st u(t 7 n)dt = u(s,n)= / [T D (l)f}(n) 



e- st g(t,l)dt 



dl 



-st 



[T D (l)f}(n)g(t,l)dl 



dt 



using ( 12. 6ft . 

By the uniqueness of the Laplace transform, 



n) 



[T D (Z)/](n)^(t,Z)dZ 



(4.11) 



f(n)h(t,X n ) 
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where h(t, A) = J °° e~ lx g(t, l)dl is the Laplace transform of E t . 

Note that inverting the </> n -transform is equal to multiplying both sides of the above 
equation by n and then summing up from n = 1 to oo. Also, the unique inverse 
u(t, .) G L 2 (D) for each fixed t > 0. Inverting the n -transform u(t,n) in (14.110 . we 
get an L 2 -convergent solution of (14.11) as 



u(t, x) 



(4.12) 



} y u(t,n)(f) n (x) 

n=l 
oo 

^f(n)(j) n (x)h(X n ,t) 



n=l 



for each t > 0. 

In order to complete the proof, it will suffice to show that the series (I4.12p converges 
pointwise, and satisfies all the conditions in (14.11) . 

Step 1. We begin showing that (14.121) convergence uniformly in t G [0, oo) in the L 2 
sense. Define the sequence of functions 

N 

(4.13) u N (t, x) = Y, f{n)(j) n {x)h{t, X n ). 

n=l 

Since g(t, I) is the density of E t , we have < h(t, X n ) < 1. Also, if < A < r/, then 
h(t,rj) < h(t, A), showing that h(t, A) is nonincreasing in A. 

Since / G L 2 (D), we can write f(x) = Yl^=i f( n )4 ) n(x), and then the Parseval 
identity yields 

oo 
n=l 

Then, given e > 0, we can choose n (e) such that 
(4.14) 



n=n (e) 



For iV > M > n (e) and t > 0, 



\u N (t,x) - u M (t,x) 



12, D 



< 



N 



E f( n ) ( l ) n(x)h(t,X r 



n=M 



2,D 



(4.15) < h(t,X no ) 2 E (f(n)) 2 < £ (/») 2 <e. 

n=n (e) n=n (e) 

Thus, the series (14.121) converges in L 2 (D), uniformly in t > 0. 

Step 2. Next we show that the initial function is defined as the L 2 limit of u(t,x) 
as t — > 0, that is, we show that t — > u(t, •) G C((0, oo); L 2 (D)) and w(t, •) takes the 
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initial datum / in the sense of L 2 (D), i.e., 

||«(V)-/lkD->0, ast^O. 

Since h(t,X) is the Laplace transform of E t , it is completely monotone and non- 
increasing in A > 0. Hence, 

oo 

u(t, x) - f(x) = f(n)(h(t, A n ) - l)0 n (x). 

71=1 

Fix e G (0, 1) and choose uq = rio(e) as in (14.14j) . Then, 

oo 

IK*, 0-/11^ = ^(/») 2 (M*,A n )-i) 2 

n=l 
no(e) 

< £(/(n)) 2 (MMn)-i) 2 

n=l 

oo 

+ E (7H) 2 (M*,a„)-i) 2 

n,=n (e)+l 

< (i-/^,A no )) 2 ||/|| 2jD + 6 

and now the claim follows, if h(t, X n ) — > 1, as t — > 0. 

This follows because £^ =>- i?o in distribution as t — > 0+ and hence the Laplace 
transforms converge. To see that E t =>- E , use the fact that {E t < x} = {W x > t} 
which is (3.16) in [32]. Then for x > and t n | 0, 

P(E tn <x) = P(W X >t n ) = l- P(W X <t n )^l- P(W X < 0) = 1, 

since W x has a density. As E t also has a density on [0, oo) we see that P(E tn < x) — 0, 
if x < 0. Thus, 

W„ < z) P(£ < x) = I(x > 0) 
at all continuity points of the limit. Note that i^o = almost surely, since W x > 
almost surely for all x > 0. 

The continuity of t i— > w(t, •) in L 2 (D) at every point £ G (0, oo) can be proved in a 
similar fashion. 

Step 3. A decay estimate for the solution u(t, x) is obtained as follows. Using 
Parseval's identity, the fact that A n is increasing in n, and the fact that h(t, A n ) is 
non- increasing for n > 1, we get 

IKv)lkD<MMi)imkD. 

Step 4. We next show that the series (14.121) defining u(t,x) is the classical solution 
to (14. ip . by proving its uniform and absolute convergence. We do this by showing 
that (14.131) is a Cauchy sequence in L°°(D) uniformly in t > 0. 
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Applying the Green's second identity, we see that 



A/(n) 



Af(x)(p n (x)dx = -X n f(n). 



i) 



Hence, Af = Y^=i —^n<Pn{x)f(n) is absolutely and uniformly convergent by assump- 
tion. Let e > 0. Since A n — > oo as n — > oo, there exists an n (e) so that for all x G D, 



(4.16) 



£ \f(n)\\Mx)\ < E A n |/»||0 n (x)|<e. 

n=no(e) n=no(e) 



This is possible since we have assumed that Af has an expansion which is uniformly 
and absolutely convergent in L°°(D). We will freely use the fact that the series 
defining / also converges absolutely and uniformly. 
For N > M > n (e) and t > and x G D, 



\u N (t,x) - u M (t,x)\ 



(4.17) 



N 

I ^2 (f) n (x)f(n)h(t,X n )\ 

n=M 
N 

< \M^\\f(n)\ < e, 

n=M 



since h(t, A) = E(e~ XEt ) < 1 for all t > and A > 0. 

This shows that the sequence u^{t,x) is a Cauchy sequence in L°°(D) and so has 
a limit in L°°(D). Hence, the series in ( 14. 12}) is absolutely and uniformly convergent. 
Also, it follows that u(t,x) satisfies the boundary conditions in ( 14. ip . 

Step 5. Next we show that the distributed-order fractional time derivative and the 
Laplacian A can be applied term by term in (14.121) . As h(t, A) is bounded above by 
unity, we have 



J2f(n)h{t,X n )A<f) n {x) 



n=l 



(4.18) 



^2f(n)h(t, \ n )\ n (t) n {x) 

n=l 

oc 

< 5Zl0n(x)||/(n)|A„ < OO, 

n=l 



where the last inequality follows from the fact that the eigenfunction expansion of 
Af converges absolutely and uniformly. Then the series 



J2f(n)h(t, A„)A0 n ( 



x 



n=l 



is absolutely convergent in L°°(D) uniformly in (0, oo). Since h(t, A) is an eigenfunc- 
tion of the distributed-order Caputo fractional derivative with ~B^h(t, A) = —\h(t, A) 
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(see O3.10p ). we have 

oo oo 

f(n)M^ {u) h(t, A n ) = f( n )Ht, A n ) A0 n ,(x). 

n=l n=l 

As the two series are equal term-by-term and the series on the right converges abso- 
lutely and uniformly, the series on the left converges absolutely and uniformly too. 

Now it is easy to check that the distributed-order fractional time derivative and 
Laplacian can be applied term by term in (I4.12p to give 



- Au(t,x) 



f(n) [<p n (x)® (u) h(t, An) - h(t, X n )A<f) n (x)} = 0, 



n=l 



so that the PDE in (14.11) is satisfied. Thus, we conclude that u defined by (14.121) is a 
classical (strong) solution to (14.11) . 
Further, using Lemma [3.21 we get 



du(t, x) 



Of 



< 



< 



El /»i 



n=l 



dh(t, X n ) 



Of 



\<t>n{x)\ 



k{t)Y,\ n \f{n)U n {x)\ := k(t)g{x). 



Since Af has absolutely and uniformly convergent series expansion with respect to 
{0„ : n G N}, we have g G L°°(D). 

Thus, it follows, from the results obtained above, that u G 7Ya(-Doo) H Cb(Doo)- 

Step 6. We next show that u G C l (D); this follows from the bounds in [2"Ul Theorem 
8.33] and the absolute and uniform convergence of the series defining /, namely, 



(4.19) 



in|l,a;r> < C(l + A n ) SUp \(f> n (x)\, 

D 



where C = C(d, A, A, dD) and A is the constant in the definition of uniform ellipticity 
of L in (TO) and A is the bound in ff2T4Tl . Here, 



u\k,a;D = SUp [Z? 7 lt] a ,D + SUp SUp |-D 7 w|, k = 0, 1, 2 



=o m=J 



and 



[D^uj^D = sup ■ 

x,y£D,x^y 
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DTm(x) -D^u(y)\ 
\x — y\ a 



are norms on C k ' a (D). Hence, 



\u 



{t,.)\l,a;D < Cy2\f{n)\h{t,\ n ){l + \ n )sUp\<P n {x)\ 



n=l 
oo 



< CV|/(n)|(l + A n )sup|0 n (x)| 
— : D 



n=l 
oo 



< 



C Y^sup \(f) n (x)\\f(n) \ + C V" A n sup \(p n {x)\\f(n)\ < oo. 

, D D 



n=l 



n=l 



Step 7. We obtain here the stochastic solution to fl4.ll) . by inverting the n -Laplace 
transform. After showing the absolute and uniform convergence of the series defining 
u, we can use a Fubini-Tonelli type argument to interchange order of summation 
and integration in the following, together with (14.121) and (14. 9p . to get a stochastic 
representation of the solution as 



u(t, x) 



POO 

Y> n (aO / WM{n)g{t,l)dl 

n=l J ° 

oo 



n=l 



g(t, l)dl 



T D (l)f(x)g(t,l)dl 



(4.20) = E x [f(X(E t ))I(r D (X)>E t )}. 

The last equality follows from a simple conditioning argument and using (12.71) . 



Step 8. Finally, we prove the uniqueness. Let u^i = 1,2, be two solutions of (14.11) 
with initial data Uj(0,x) = f(x) and Dirichlet boundary condition Ui(t,x) = for 
x G 3D. Then U = U\ — Ui is also a solution of (14.11) with zero initial data and zero 
boundary value. Taking n -transform on both sides of (14. ip we get 



6(") 



U(t, n) = -X n U(t, n), U(0, n) = 0, 



and then U(t,n) = for all t > and all n > 1. This implies that U(t,x) = in 
the sense of L 2 functions, since {0 n : n > 1} forms a complete orthonormal basis for 
L 2 (D). Hence, U(t, x) — for all t > and almost all x G D. Since U is a continuous 
function on D, we have U(t,x) = for all (t,x) G [0, oo) x D, thereby proving the 
uniqueness. □ 

Corollary 4.2. The solution in Theorem ( f^.j| ) also has the following representation: 

u(t,x) = E x [f(X(E t ))I(T D (X)>E t )]=E x [f(X(E t ))I(T D (X(E))>t)}. 
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Proof. The argument is similar to [3"4"l Corollary 3.2] and so we only sketch the proof. 
Given a continuous stochastic process X t on M. d , and an interval / C [0, oo), we denote 
X(I) = {X u : u G I}. Since the domain D is open and X u is continuous, it follows 
that 

{r D (X)>t} = {X([0,t])cD}. 

Next note that, since E t is continuous and monotone nondecreasing, E([Q, t}) = [0, E t }. 
Finally, we observe that 

{t d (X(E)) > t} = {X(E([0,t}) C D} 

= {X([0, E t ]) CD} = {r D (X) > E t } 

which completes the proof. □ 

Remark 4.3. If we time-change Brownian motion B(t) using a nondecreasing stable 
Levy process Wt, then the conclusions of Corollary 14.21 do not hold, since the stable 
subordinator does not have continuous sample paths; see, for example, Song and 
Vondracek [16]. In our case, killing the process B(t) and then applying the time 
change E t is the same as applying the time change and then killing, since the sample 
paths of E t are continuous. 

The next result establishes existence of strong solutions of distributed-order frac- 
tional Cauchy problems (13.71) with L = A. 

Corollary 4.4. Let f e C^ k (D) be a Ik-times continuously differentiable function of 
compact support in D. If k > 1 + 3d/ 4, then &4.1\ ) has a classical (strong) solution. 
In particular, if f G C^°(D), then the solution of d4-l\) is in C°°(D). 

Proof. By Example 2.1.8 of [16], \<j) n (x)\ < (\ n ) d/4 - Also, from Corollary 6.2.2 of [17], 
we have \ n ~ n 2 l d . 

Applying Green's second identity k-times, we get 



(4.21) 



A*/(n) 



k k f(x)<p n (x)dx = (-\ n ) k f(n). 



i) 



Using Cauchy-Schwartz inequality and the fact / G C\ (D), we get 



A*/(n) < 



(A k f(x)) 2 dx 



D 



1/2 r 



■n 1/2 



(4> n (x)) 2 dx 



D 



(f 2k (x)Ydx 



D 



1/2 



where is a constant independent of n. 
This and fl4T2T|) give \ f(n)\ < c k (\ n )- k . 
Since 



A/(x) = ^-A n /»0 n i 



x) 



n=l 
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to get the absolute and uniform convergence of the series defining A/, we consider 

oo oo 

^A n |0 n (x)||/»| < J2( X n) d/i+1 C k (X n )- k 
n=l n=l 

oo oo 

< c k Y,(n 2/d ) d/i+1 - k = c k J2^ 1/2+2/d ~ 2k/d 

n=l n=l 

which is finite if - | - \) > 1, i.e., if k > 1 + \d. □ 

Remark 4.5. In an interval (0, M) C M, eigenf unctions and eigenvalues are explicitly 
known. Eigenvalues of the Laplacian on (0, M) are (ntr/M) 2 , and the corresponding 
eigenfunctions are sin(n7rx/M), for n = 1, 2, ■ • • . The form of the solution in ( 14.121) 
on a bounded interval (0, M) in K. was obtained by [TJ [37] . Agrawal pQ worked with 
single-order fractional Cauchy problem with Dirichlet boundary conditions. Naber 
considered Dirichlet and Neumann boundary conditions in one space dimension. 



Recall that = (0, oo) x D and define now 

HlODoo) = {u:D 00 ^R: Lu{t, x) G Cpoo)}; 

H b L (D 00 )=n L (D 00 )n{u:\d t u(t,x)\<k(t)g(x), geL°°(D), t > 0} 

where k(t) is defined in ( 13.171) . The following result extends Theorem 14. II to general 
uniformly elliptic second-order operators. 

Theorem 4.6. Let D be a bounded domain with 3D e C 1,a for some < a < 1, 
and suppose that L is given by ( 12.11) with aij G C a (D). Let {X(t)} be a continuous 
Markov process with generator L, and To{t) the killed semigroup corresponding to the 
process {X(t)} in D. Let E t be the inverse (13. 3p of the subordinator W t , independent 
of {X(t)}, with Levy measure (13. 5p . Suppose that fi(d/3) = p(j3)dj3, as in Lemma 
\3.2\ and is the distributed- order fractional derivative defined by (13.81) . Then, for 
any f G D(Ld) PI C l (D) PI C 2 (D) for which the eigenfunction expansion (of Lf ) with 
respect to the complete orthonormal basis {ip n : n > 1} converges uniformly and 
absolutely, the (classical) solution of 

(4.22) B {u) u(t,x) = Lu(t,x), x G D, t > 0; 

u(t,x) = 0, x G 3D, t > 0; 
u(0, x) = f(x), x G D, 
for u G HKDoo) n C 6 (.Doo) n C\D), is given by 

u(t,x) = E x [f{X{E t ))I{r D {X) > Et)] = E x [f{X{E t ))I{r D {X{E)) > t)} 

(4.23) = / T D {l)f{x)g{t ) l)dl = y^f{n)^ n {x)h{t^ n ) ) 

Jo o 

where h(t, /i) = E(e -MjEt ) = J °° e~ Xfl g(t, x) dx is the Laplace transform of E t . 
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Proof. Let u(t, x) = G(t)F(x) be a solution of fTP2|) . Substituting in the PDE fTQ2D 
leads to 

F(x)D M G(t) = G(t)LF(x) 
and now dividing both sides by G(t)F(x), we obtain 

B^G(t) _ LF(x) 
~G{t)~ ~ ~F(x)~ ~ ~^ 

That is, 

(4.24) B (u) G(t) = -fiG{t), t > 0; 

(4.25) LF(x) = -fiF(x), x e D, F\ 9D = 0. 

Problem (14.251) is solved by an infinite sequence of pairs (fi n ,ip n ), n > 1, where 
< fii < /i 2 < fa < ■ ■ ■ is a sequence of numbers such that fi n — > oo, as n — > oo, 
and ip n is a sequence of functions that form a complete orthonormal set in L 2 (D) (cf. 
(12.51) ). In particular, the initial function / regarded as an element of L 2 (D) can be 
represented as 

oo 

(4.26) /(z)=$^/(nM>(z). 

n=l 

An application of the Parseval identity yields 

oo 

(4-27) ll/l| 2 ,D = ^(/(n)) 2 . 

n=l 

Using the \x n determined by (I4.25P and recalling from Lemma I'd. II that h(t,/j, n ) solves 
(I4.24p with n = fi n , we obtain 

G(t) = G (n)h(t lf i n ), 
where Go(n) is selected to satisfy the initial condition /. We will show that 

oo 

(4.28) u(t, x) = J2 f(n)h(t, ^ n {x) 

n=l 

solves the PDE (T4T221 . 

Define approximate solutions of the form 

N 

(4.29) u N {t, x) = Y, G Q (n)h{t, fi n )Mx), G {n) = f{n). 

n=l 

Step 1. Following the proof of Theorem 14.11 it can be shown that the sequence 
{u]y(t, -)} N< z N is a Cauchy sequence in L 2 (D) PI L°°(D), uniformly in i 6 [0, oo). 
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Hence, the solution to (14.221) is given formally by 

oo 

(4.30) u(t, x) = Y, f(n)h(t, Hn)if> n (x)- 

n=l 

Step 2. It follows again, as in the proof of Theorem 14. 1[ the series defining u and Lu 
converge absolutely and uniformly so that we can apply the fractional time derivative 
and uniformly elliptic operator L term by term to show that u defined by (I4.30p is 
indeed a classical solution to (14.221) . 

Step 3. The stochastic representation of the solution, as given in (I4.23p . also follows 
in a similar manner and hence we omit the details. 

Step 4. We have also a decay estimate for u as in Theorem 14. 1[ namely, 

IK-,t)|| 2 ,D<M^i)ll/lkz> 

The uniqueness of the solution can be proved as before. □ 

5. Extensions and open questions 

In this section, we derive some general conditions on the mixing distribution fi(dj3) 
in (13.81) that are sufficient to obtain classical solutions to the distributed-order frac- 
tional Cauchy problem (13. 7p on bounded domains, as in Theorems 14.11 and 14.61 In 
particular, we remove the assumption that n(d[3) = p{(3)d(3 to allow atoms. Then we 
discuss related literature, and some open questions. 

Lemma 5.1. Suppose fi is a finite measure with supp(fi) C (0, 1) that satisfies (13.91) . 
Assume also that \dth(t, A)| < b(\)k e (t) for some functions b and k e satisfying the 
condition 

(5.1) 6(A)^ 1 ^^|^)<oo, 

for t, A > 0. Then h(t, A) = K(e~ XEt ) is a classical solution of the eigenvalue problem 



(5.2) B^h{t, A) = -Xh(t, A); h(0, A) = 1. 

Proof. The proof follows from Lemma 13.11 and the fact that ( 15.11) is a sufficient 
condition for D^/i(t, A) to be defined as a classical function. □ 

Suppose that \d t h(t, A)| < b(\)k e (t), where b and k e satisfy, in addition to ( 15. II) . 

oo 

(5.3) Mt)5>(A n )/(n)|0 n (x)| < oo. 

n=l 

It is assumed here that the above series converges absolutely and uniformly for t > 0. 
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Let (3 £ (0, 1), -Doo = (0, oo) x D and define 

Wi(Doo) = {u : R : s) £ C(A*)}; 

H b £%D x )=H L {D 0O )n{u:\d t u(t,x)\<k e (t)g(x), g £ L°°(D), t>0}, 

where k e and 6 satisfy (15.11) and (15.31) . The following result extends Theorem 14.61 to 
allow atoms in the mixing measure fi(d/3). 

Theorem 5.2. Let D be a bounded domain with 3D £ C 1,a for some < a < 1, 
and suppose that L is given by (12.11) with aij £ C a (D). Let {X(t)} be a continuous 
Markov process with generator L, and To(t) be the killed semigroup corresponding 
to the process {X(t)} in D. Let E t be the inverse (13. 3p of the subordinator Wt, 
independent of{X(t)}, with Levy measure (1331) . Let f £ D(L D )nC 1 (D)nC 2 (D) for 
which the eigenfunction expansion ( of Lf ) with respect to the complete orthonormal 
basis {ip n : n > 1} converges uniformly and absolutely. Then the (classical) solution 
of 

(5.4) B {u) u(t,x) = Lu(t,x), x £ D, t > 0; 

u(t,x) = 0, x £ 3D, t > 0; 
u(0,x) = f(x), x E D, 

for u £ H^ e (D 00 ) fl Cb^Doo) fl C l (D), with the distributed order fractional derivative 
defined by (13.81) . zs gwen fry 

u(t,a;) = E x [f(X(E t ))I(m(X) > E t )} = E x [f(X(E t ))I(r D (X(E)) > t)} 

POO 00 

(5.5) = / T D {l)f{x)g(t,l)dl = Y^f{n)^ n {x)h(t,ii n ), 

Jo o 

where h(t,/i) = K(e~^ Et ) = J °° e~ Xfl g(t, x) dx is the Laplace transform of E t . 

Proof. The proof follows the same steps as in the proof of Theorems 14.11 and 14.61 and 
using the properties (15.11) and (I5.3P of the function k e (t). □ 

Remark 5.3. Let fi(d/3) = E^Li^W 1 ~ Pj))~% s (P) d P> for < ft < ft < ■■• < 

j3j\j < 1. In this case, the subordinator is W t = YljLi c 3^t 1 f° r independent stable 

subordinators W t \ for j — 1, • • ■ , N. In this case, the functions k e (t) and 6(A) that 
satisfy 1ET3J) . flO) and 

|5 t /i(t,A)| <b(X)k e (t) 

are fc e (t) = (c^ J sin(/5j7r)) _1 (t^ _1 ) for all j = 1, ■ ■ • , iV and 6(A) = A, respectively. 
The proof of this fact follows the same steps as in the proof of equation (2.19) in [27] 
using the properties of /x(/3). Hence, in this case Theorem 15.21 applies and we have a 
classical solution of (15.41) given by (15. 5p . 
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To conclude this paper, we discuss the related literature and some open problems. 
Given a uniformly bounded, strongly continuous semigroup T(t) with generator L on 
some Banach space H, the Cauchy problem du(t, x)/dt = Lu(t, x) with u(0, x) = f(x) 
has solution u(t,x) = T(t)f(x) for any initial condition / e H [3\. One important 
special case is the pseudo-differential operator 



that appears in the backward equation of a Markov process X(t) [2~4"1 |4"4"] . The 
probability distribution of the Markov process X(t) solves the forward equation, which 
is the Cauchy problem with the adjoint of the generator L. The integral term in (15.61) 
represents a jump diffusion (e.g., a stable process). For stable generators, the explicit 
connection with stochastic differential equations driven by a stable Levy process was 
established by Zhang et al. [H] and Chakraborty In that case, the integral term 
in (15.61) can be written in terms of fractional derivatives in the space variable. 

The solution (14. 2 j) to the fractional Cauchy problem d^u(t,x)/dt^ = Lu(t,x) for 
< (3 < 1 was established by Baeumer and Meerschaert [I] in the general Banach 
space setting. Baeumer et al. [5] and Nane [3S] specialized to Markov processes with 
generator (15.61) . and established a connection to iterated Brownian motion [21 El El 
ITS] . Hahn et al. [22] developed the connection with stochastic differential equations 
driven by a time-changed Levy process X(E t ) for generators (15.61) . so that their result 
includes jump diffusions on M. d . Their results extend those of [TTl HSJ to distributed- 
order fractional Cauchy problems on M d . Hahn et al. [22] also give the integral solution 
(14.21) as in [HE]. Kochubei [27] provides strong solutions of distributed order fractional 
Cauchy problems on M. d in the case L — A. Meerschaert and Schemer [33] discuss 
generalized Cauchy problems of the form ijj w (d/dt)u(t, x) = Lu(t, x) + 5(x)ipw{'ti °°), 
where ipw{ s ) is the Laplace exponent of a nondecreasing Levy process (subordinator) 
whose Levy measure (pw has infinite total mass, and L is the generator of another Levy 
process. This reduces to the distributed order fractional Cauchy problem (13.71) in the 
special case when (13.41) holds. As in Section [31 the paper [33] shows that the density 
u(t, x) of the CTRW scaling limit A(E t ) solves the generalized Cauchy problem, when 
E t is the inverse of the subordinator W t with E[e _sH/ *] = e _i ^ w ^- ) . Strong solutions of 
generalized Cauchy problems on IR d , including fractional or distributed-order Cauchy 
problems, seem to be an open problem. 

For bounded domains, the general results of [1] remain valid, so that the solution 
formula (14.21) still holds in the appropriate Banach space. The results of this paper 
provide strong solutions in that case, so long as L generates a diffusion without 
jumps. To the best of our knowledge, the construction of strong solutions for jump 
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(5.6) 




diffusions remains a challenging open problem. Eigenvalue expansions can be found 
explicitly in some special cases. The main technical difficulty is to obtain regularity 
of the eigenfunctions, or at least sharp bounds, for the generator (15.61) in the case of 
jump diffusions on bounded domains. See Chen et al. [H] for a recent study on this 
problem. One explicit example is to take L = — (— A) Q//2 for < a < 2, the classical 
fractional power of the Laplacian [23], which generates a spherically symmetric stable 
Levy process. This results from (15.61) with a = b = and (f)(x,dy) = Cd, a \\y\\~ a ~ 1 dy, 
where Cd, a is a constant that depends on the stable index a and the dimension d of 
the space, see for example [29] . This is a type of fractional derivative in space, called 
the Riesz fractional derivative of order a. Some results for this case are available in 
Chen and Song [12], Chen et al. [13] and Song and Vondracek [16]. 
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